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Abstract In this paper for studying the anisotropic strange quark stars, we 
assume that the radial pressure inside the anisotropic star is a superposition of 
pressure in an isotropic case plus a Gaussian perturbation term. Considering a 
proportionality between electric charge density and the density of matter, we 
solve the TOV equation for different cases numerically. Our results indicate 
that anisotropy increases the maximum mass Mmax and also its corresponding 
radius R for a typical strange quark star. According to our calculations, an 
anisotropy amplitude of A = 3 x with a standard deviation of 

cr = 3 X lO^m leads to a neutron star of Furthermore, electric charge 

not only increases the maximum mass and its corresponding radius, but also 
raises up the anisotropy factor. We can see that the tangential pressure pt 
and anisotropy factor A unlike the radial pressure Pr have a maximum on the 
surface and this maximum increases by adding electric charge effect. However, 
we show that anisotropy can be more effective than electric charge in rasing 
maximum mass of strange quark stars. 
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1 Introduction 

Strange quark stars are hypothetical type of compact exotic stars which are 
firstly speculated in Refs. Describing the structure of stars, including 

mass (M) and radius (i?) will be available by solving hydrostatic equilib¬ 
rium equations. In the case of compact stars, due to high density. General 
Relativity (GR) dominates and Newtonian hydrostatic equilibrium should 
be replaced by its GR counterpart. Assuming a static and spherical sym¬ 
metric geometry and an isotropic matter, Einstein field equations lead to 
Tolman—Oppenheimer—Volkoff equation (TOV) [1]. The structure of com¬ 
pact stars has been investigated by Refs. [siimmis] and many others by solving 
TOV equation using a suitable equation of state. Now in this paper, we deal 
with the strange quark matter (SQM) which its EOS is well described by MIT 
Bag Model [S]. 

nuclear matter in very high densities can be anisotropic and it is a great mo¬ 
tivation to study the effect of anisotropy in the structure of relativistic stars 
[SIS]. Additionally, existing a solid core and a strong magnetic field in neu¬ 
tron stars can be related to anisotropy in the matter of star pfTO] . Solving 
anisotropic TOV equation requires a physically reasonable assumption. Eor in¬ 
stance, in Ref. [B] a specific density profile p{r) has been chosen and in Ref. [5] 
a special metric function A{r) has been utilized to solve it analytically. But 
in this paper, we try to apply a numerical method for solving modified TOV 
equation by considering a perturbation term. Initially, we consider an isotropic 
and free charge matter and obtain a solution for radial pressure Prij")- Then 
in anisotropic case we use the previous Pr(f’) as the unperturbed solution in 
order to study R and M in this kind of stars. It should be noted that in this 
work we are dealing with a non-rotating strange quark star but we can expect 
that rotation allow a larger maximum mass (about 40%) according to Ref. [5] . 
This paper is organized as follows: In §2 we first solve the TOV equation for 
an isotropic matter and then compare it with a charged and anisotropic one. 
Discussion and concluding remarks are represented in §3. 


2 TOV Equation 


For describing mass and radius of a self gravitating configuration in relativistic 
term we have to consider a suitable metric. The line element in interior of star, 
assuming a static and spherically symmetric geometry can be written as: 

ds^ =-c^e^^dt^+ e^'^dr^+r^dn^, ( 1 ) 


where c is the light speed while (j) and A are spherically symmetric metric 
functions. 

In order to obtain a General Relativistic hydrostatic equilibrium, we have to 
solve the Einstein equations: 


= 


SttG 


_ T 

4 


( 2 ) 
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for the metric given in Eq. ©• Also we must choose a reasonable energy- 
momentum tensor which satisfies the following conservation law: 

= 0. (3) 

Now in the next subsections we try to obtain the TOV equation for two dif¬ 
ferent cases. 


2.1 Uncharged Isotropic Matter 


First of all we consider an uncharged perfect fluid. The energy-momentum 
tensor of a perfect fluid is given by m-- 


Tlf = 


f-pc? 0 0 0\ 
0 pOO 
0 Op 0 
V 0 oopj 


( 4 ) 


where p and p are the pressure and mass density respectively which are 
spherical symmetric [5] . We also define the gravitational mass as [T^ : 


m{r) = / p{x)dx. (5) 

Jo 

By using the above equation, the Einstein equations ([5]) and the conservation 
law for metric o yield to hydrostatic equilibrium equation: 


dp 

dr 


2 4:'!TGpr^ + uiGc^ 

rc2(rc2 -2mG) ' 


( 6 ) 


Noticeably m, p, and p are functions of radial coordinate (r). Equation ([S|) 
together with Eq. ([B]) is called TOV equation which for its solvability it has 
to be supplemented by an EOS m- In this work we consider the MIT Bag 
Model for EOS of SQM as: 


^ = a{p-ps). (7) 

According to Ref. [2] for SQM2 model ps = 3.056 x 10^'^kg m~^ and a = 
0.324 . Now we can solve TOV equation numerically, using physical boundary 
conditions m- 


p(0) = Pc 

(8a) 

m(0) = 0 

(8b) 

piR) = 0 

(8c) 

m{R) = M 

(8d) 


In addition, in the stellar interior, following conditions should be satisfied: 
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p > 0 

(9a) 

dr 

(9b) 

dp 

(9c) 

±>0 

dp 

(9d) 


Conditions (li^ and ([90 are trivial for preserving hydrostatic equilibrium 
and (1^ refers to causality inside the star which does not permit a > 1 in Eq. 
© 0. Also the EOS should satisfy the microscopic stability condition (IMl) . 
otherwise the star collapses spontaneously [T5] . In the next section we will use 
the obtained pressure p{r) as the unperturbed radial pressure for anisotropic 
case. 


2.2 Charged Anisotropic Matter 


In this section instead of a perfect fluid we are dealing with an anisotropic 
matter. It means that pressure in radial and tangential directions are not 
necessarily the same. With this assumption, which is predicted in very high 
density ranges0, the energy momentum tensor in CCS unit system will be |51 

HE]: 


^-(pc2 + ^) 0^0 0 \ 

0 Pr-^ 0 ^ 0 

0 0 Pt + ^ 0 

V 0 0 ^ Pt + ^J 


( 10 ) 


where pt and Pr are tangential and radial pressure respectively. The gravita¬ 
tional mass which is the total contribution of the energy density takes the new 
form as m- 

m(r) = J Airx'^ (^p{x) + dx, (11) 


where E is the radial electric field, defined as: 


E{r) = \ [ 'inx'^pche'^dx (12) 

Jo 

and pch is the charge density. Defining A = pt — Pr as anisotropy factor 0 , 
Einstein equations (0) and conservation equation (0) lead to: 


e 


-2A 


= 1 - 


2mG 
rc^ ’ 


(13) 
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dpr 

dr 



4:TTGr^{pr — ■^) + rnGc^ 
rc^{rc^ — 2mG) 


+ Epche^. 


(14) 


It should be noted that if Z\ > 0, in Eq. (O the magnitude of radial pressure 
derivative decreases. This causes a softer falling of radial pressure which results 
higher maximum masses and their corresponding radii. Equations (ED and 
(fTTl) represent the modified form of TOV equation. It is easy to see that putting 
E = 0 and Z\ = 0 in Eq. ED, simply return uncharged and isotropic cases. In 
an anisotropic matter, MIT Bag Model is [6]: 


= a{p- ps). 


(15) 


We should notice that tangential pressure pt does not necessarily vanish on the 
surface of star, in contrary to Pr.[6j|8]. Solving modified TOV equation requires 
some extra assumptions. Assuming that anisotropy has a slight effect on 
we add a perturbation to isotropic pressure profile to obtain pr as: 


Prip^Anisotropic — p{‘^)Isotropic E Spip'). (ff^) 

The following perturbation function has to satisfy these boundary conditions: 


Pr(0) =Pc, 

(17a) 

Pr{R) = 0, 

(17b) 

— U, 
dr 

(17c) 

dPr\ „ 

—^\r=R = U. 

dr 

(17d) 


Regarding to Eq. [HI p{r)isotropic and Pr{r) Anisotropic have to satisfy men¬ 
tioned boundary conditions. So Sp{r) must meet those conditions too. Gaus¬ 
sian function is the simplest choice (but not the only on^il) for fulfilment of 
the conditions I17al through Il7dl by setting the position of Gaussian pressure 
prohle at /i = ^ in the following equation: 

dp(r) = Aexp . (18) 

One can see that this Gaussian perturbation can satisfy mentioned boundary 
conditions physically. We need 6 and its derivative to be nearly zero at origin 
and also on the surface. In fact we can tune A and a so that the perturbation 
to be so small at those points. According to Ref. the anisotropy can cause 
an increase in the radial pressure. Hence we have chosen a positive and small 
enough value compared to Pr for parameter A. Eventually studying charged 


^ For example Lorentzian function could be another choice for Sp{r) which yields no 
significant difference. 
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case requires one more physically assumption. According to Ref. electrical 
charge density can be related to the density of matter as follows: 


Pch = f X P, (19) 

which is a physically reasonable assumption[T7]. In this equation / is the 
charge fraction that can be considered / < 10“^esM g~^ to satisfy causality, 
stability, and electrical neutrality conditions of stars [12]. Therefore in the 
most general case (charged anisotropic matter), there is 7 variables namely 
m,pr,pt, p, Pch, E, A and 7 equations (ITTl) - (fTOl) and (|TO)) . We have used Runge 
Kutta Fehlberg fourth-fifth order method (RKF45 Method) to solve the cor¬ 
responding coupled ODE systems in each cases below: 

1. Isotropic uncharged (Z\ = 0, A = 0 ), 

2. Isotropic charged {A = 0, E ^ 0), 

3. Anisotropic uncharged {A ^ 0, E = 0), 

4. Anisotropic charged {A ^ 0,E ^ 0). 


3 Discussion 

In this paper unlike the analytical solutions for anisotropic matter such as 
Ref. [5] and Ref. jS] we have used numerical result of isotropic solution and a 
Gaussian perturbation rather than using a totally assumed function. In the fol¬ 
lowing subsections we will discuss the result of adding anisotropy and electric 
charge to TOV equations. Our solution indicate that adding positive Gaussian 
anisotropy and electric charge to TOV equation increase the maximum mass 
of neutron star and can predict more massive neutron stars in agreement with 
Ref. [TH| . 

3.1 Gharge Effect 

Solving modified TOV equation reveals that adding electric charge makes the 
radial pressure Pr to be increased, as we can see in Eig. [TJ But the sensitivity 
of tangential pressure pt to the electric charge is much more, so A raises up in 
charged case as is apparent from Eig. [2| This behavior is not only the direct 
result of our speculation about Pr, but also comes from the nature of TOV 
equation. 

Fig. [3] and Fig. H] refer to uncharged and charged isotropic cases. One may 
expect adding electric charge to TOV equation causes an increase in the max¬ 
imum mass of strange quark star Mmax and also in the corresponding radius 
R for obtained Mmax- Since in GR gravitational mass of star and total en¬ 
ergy are proportional, adding charge density pcu increases the total energy 
and therefore the above result makes sense. [T3]. However, electrical neutrality 
condition of stars does not allow a significant growth in maximum masses and 
their corresponding radii. 
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Fig. 1 Radial pressure pr for various cases: isotropic uncharged, isotropic charged, and 
anisotropic uncharged matter with / = 5 x 10“^esu g~^, A = 3x cr = 3x lO^m 

and p = ^ versus radial coordinate r, all having the same central density of pc = 1 x lO^^kg 
m-3. 


3.2 Anisotropy Effect 

It is shown in Fig. 0 and Fig. [5] that the maximum mass of strange quark star 
has been raised by adding Gaussian perturbation for obtaining pr. Further¬ 
more, according to Fig. [5] the tangential pressure pt and anisotropy A have a 
maximum and they do not vanish on the surface. In fact our solution indicate 
that if there is an anisotropy in the star, it should be maximum on the surface. 
Recent measurements indicate that there exist a pulsar of mass 1.97 ± 0.04 Mq 
and this mass rules out nearly all currently proposed equations of state [18]. 
We have used a trial and error method to get the mentioned mass above. Our 
calculations indicate that an anisotropy amplitude of A = 3 x lO^^iVm”^ with 
a standard deviation of cr = 3 x lO^m and p = ^ in Eq. (fT51) can survive SQM 
equation of state which satisfies boundary and hydrostatic equilibrium condi¬ 
tions by 1% and 10% uncertainties at the origin and on the surface respectively 
as is depicted in Fig. |T] 


4 Summary 

We have studied anisotropic charged strange quark stars with a Gaussian 
perturbation. Our calculations have shown that the electrical charge effect is 
much less effective in increasing the maximum mass of star in comparison to 
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Fig. 2 Anisotropy factor A for anisotropic charged and uncharged matter with / = 5 x 
10“^esu g~^ versus radial coordinate r with a same central density pc = 1 x lO^^kg m~^. 



Fig. 3 Gravitational mass M versus central density pc in two cases: Uncharged (solid curve) 
and charged isotropic (dashed curve) matter. 
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Fig. 4 Gravitational mass M versus radius of star R in two cases: Uncharged isotropic (solid 
curve) and charged isotropic (dashed curve) matter. 



Fig. 5 Gravitational mass M vs central density pc for uncharged isotropic (solid curve) and 
uncharged anisotropic (dashed curve) matter. 
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Fig. 6 Gravitational mass M vs radius of star R in two cases: Uncharged isotropic (solid 
curve) and uncharged anisotropic (dashed curve) matter. 


anisotropy effect. Also our work has supported that anisotropy can be one of 
the candidates of describing massive neutron stars with current equations of 
state. 
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